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Abstract. We probe the region around the protostar 
HH108MMS by deep mid infrared photometric and po- 
larimetric imaging. The protostar is detected at 14/im 
in absorption against the diffuse background. Next to 
HH108MMS, we find a second absorbing core, named 
Ql, and the young stellar object IRAS18331-0035 which 
is more advanced in its evolution and already seen 
in emission at 12/im and 14/im. HH108MMS, Ql and 
IRAS 1833 1-0035 form a triplet along an extended filamen- 
tary absorption feature. From the variation of the surface 
brightness across the source, we derive for HH108MMS 
and Ql the optical depth and density profile. Along the 
axes which are parallel to the filament, the density distri- 



butions follow a p oc 



^-1.8 



power law. We estimate that 



the intensity of the background radiation at 14/im is about 
two times stronger than the intensity of the interstellar ra- 
diation field in the solar neighborhood. The present photo- 
metric data of IRAS18331-0035 between 12/im and 1.3mm 
can be explained by a central source with a luminosity of 
2.5 L0 that is surrounded by a spherical cloud of 1.1 M© 
with a 1/r density distribution. As HH108MMS is also 
seen in the millimeter dust emission, we can derive the ra- 
tio of the dust extinction coefficients at 14/im and 1.3mm 
and obtain Ki4^m/Ki3oo;im ~ 470. Because models for the 
dust in the diffuse interstellar medium predict a ratio of 
around 2000, our value points to fluffy composite grains 
which are expected to prevail in dense and cold environ- 
ments. 

First mid infrared polarisation images of pre-stellar ab- 
sorbing cores are presented. At 12/im and 14/xm the polari- 
sation of the region around HH108MMS is strong (> 15%) 
and tightly correlated with the source triplet. We demon- 
strate that the high degree of polarisation can be explained 
by extinction of rotationally aligned dust particles of mod- 
erate elongation. 



Send offprint requests to: rsiebenni@eso.org 

* Based on observations with ISO, an ESA project with in- 
struments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United 
Kingdom) with the participation of ISAS and NASA. 



Key words: ISM: dust, extinction - Stars: pre-main se- 
quence - Stars: individual: HH108MMS - Techniques: po- 
larimetric - Infrared: ISM: continuum 



1. Introduction 

Herbig Haro objects mark shock regions in the outflow 
from very young stars. The 1.3mm continuum survey of 
Herbig Haro objects by Reipurth et al. (1993) indicates 
that the young stars, which are the energy sources of the 
outflow, are surrounded by dusty envelopes of about 0.1 
- 3 solar masses. The circumstellar material is so massive 
that the stars are probably still in the accretion phase. In 
a follow-up study of the 1.3mm dust continuum, Chini et 
al. (1997) discovered in the Serpens star forming region 
which is at a distance D — 310 pc (De Lara et al. 1991) 
a protostellar candidate located 71" (0.1 Ipc) away from 
the source IRAS18331-0035. The latter is beheved to be 
the driving engine of the two aligned nearby bow shocks 
HH108 and HH109 (Reipurth & Eiora 1992, Ziener & 
Eisloffel 1999); they he to the South- West of IRAS18331- 
0035, about O.lpc and 0.21pc away. In this paper, we 
present photometric and polarimetric mid infrared images 
of the region. 

2. Observations 

2.1 Photometric imaging 

The protostellar system HH108MMS was photometrically 
imaged on 5 March 1996 with the long wavelength array 
of ISOCAM (Cesarsky et al. 1996) on board the infrared 
space observatory (ISO). The morphology of the object 
was studied by applying the observing template CAMOl 
in raster mode. A 5 x 5 micro-raster was performed with 
step size of 20" in both space craft directions; we used 
the broad band filter lw3 (12 - 18 nm), the 6" lens and 
the large field mirror. After 100 stabilization frames, 33 
exposures were read out at each raster position. The read- 
out time for each exposure was 2.1s. 

The data were reduced with the ISOCAM interactive anal- 
ysis (CIA). The basic reduction steps such as dark current 
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Fig. 1. ISOCAM lw3 image (14.3^ni) in grey scale around the protostar HH108MMS (30" 9300AU). The contour 
Hnes are an overlay from Chini at al. (1997) and show the 1.3mm dust emission. In the mid infrared, the right 
1.3mm source, IRAS18331-0035, is in emission, the left, HH108MMS, is seen in absorption. The other absorbing 
knot in between is denoted Ql; it is weak at 1.3mm, but clearly detected. There is a hint of a third absorbing core 
in the uppermost West. The ISOCAM image shows a large filamentary structure into which IIH108MMS, Ql and 
IRAS18331-0035 are embedded. 



subtraction, initial removal of cosmic ray hits, detector 
transient fitting and photometry are described in Sieben- 
morgen et al. (1999). Because in our observing strategy 
several detector elements saw the same part of the sky 
at different times, we used the redundant information to 
correct for the long term camera drift and to fine-tune 
the deglitching process (Miville-Deschenes et al. 2000). 
The coadded images at each raster position were then 
sky-projected and corrected for field distortion. The fi- 
nal raster map has a pixel scale of 2" which corresponds 
to ^620AU at the distance of Serpens {D = 310 pc). 
Fig. |l| displays in grey scale our lw3 image of HH108MMS 
and its surroundings. Superimposed are contours of the 
1.3mm emission as observed by Chini et al. (1997). The 
brighter 1.3mm source coincides with IRAS 18331-0035. 
At 14/xm we detect this source in emission, HH108MMS 
itself is seen in absorption. Besides a second absorption 
core, Ql, and a tentative third absorbing knot, we find 



a filamentary structure again in absorption against the 
diffuse infrared background. The filament extends about 
200" (0.3 pc) along its major and 20" (0.03 pc) along its 
minor axis. The major axis is oriented at a position angle 
PA = 45°. 

Photometry of IRAS18331-0035 in a 12" aperture gives 
Fi4f_i„i = 4.2 ± 0.3 mJy. The fluxes have been scaled by a 
correction factor for a monochromatic point spread func- 
tion and color-corrected assuming a i^A oc A spectrum in 
the band. 

2.2 Polarimetric imaging 

Polarimetric images of the HH108MMS region were ob- 
tained on 30 March 1998 with ISOCAM. The polarisation 
pattern was studied by applying the observing template 
CAM05 (Siebenmorgen 1996) in broad band filters lw3 
(14.3 /im) and IwlO (12.0 ^m). In each filter, a 2 x 2 
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Fig. 2. Polarisation map at 14.3/ini (lw3) of the protostel- 
lar system HH108MMS (cmp. with Fig. 1). The intensities 
are shown as contour overlay and polarisation vectors as 
bars. For the absorbing cores the polarisation vectors in- 
dicate the magnetic field direction. 
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Fig. 3. As Fig. | for filter IwfO (12.0 ^im). 



micro-raster was performed at a step size of 54" in both 
space craft directions. The camera was configured in the 
6" lens and the large field mirror. After 30 stabilization 
frames, 15 exposures were read out at each raster posi- 
tion. The read-out time for each exposure was 5s at gain 
two. After a raster without polarisers, similar rasters were 
taken for each of the three polarisers. The polariser rasters 



1.40 
1.30 ^ 

1.20 



« I.IOH 

P : 



1.00 

0.90 



P3 



Ict 




Fl 



cycle 

Fig. 4. Flux of IRAS18331-0035 as measured through fil- 
ter lw3 (14.3 /im) and polariser PI (*), polariser P2 {()) 
and polariser P3 (A) for the different observing cycles. 
The typical Icr error bar is indicated. 
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Fig. 5. As Fig. | for filter IwlO (12.0 /im). 

were repeated over three observing cycles. The nominal 
pointing position (18*^ 35'° 46".5 - 0° 32' 51", J2000.0) of 
the target was corrected in the polariser rasters for the 
known source displacement on the detector. 
The polarisation maps of HH108MMS in filters lw3 and 
IwlO are shown in Fig. ^ and Fig. ^. Superimposed are con- 
tours of the intensities without the polarisers. They reveal 
the same structure as detected in the CAMOl observa- 
tions. Photometry of IRAS18331-0035 at 12 fim (IwlO) 
gives 2.0 ± 0.2 mJy and at 14 fiui (lw3) 4.2 ± 0.4 mJy, in 
accordance with the CAMOl observations. 
As in the CAMOl observations described above, we ap- 
plied the basic data reduction steps, such as dark current 
subtraction, initial removal of cosmic ray hits and tran- 
sient correction. A flat field is derived by taking advan- 
tage of the fact that all detector pixels spend more time 
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Fig. 6. A sketch showing how we envisage that 
HH108MMS is embedded into the molecular cloud. We 
indicate the temperature and optical depth in the various 
regions. Looking towards HH108MMS, one observes the 
intensity Is, away from it /b; Iq and /z are the intensity 
behind the Serperns cloud and of the foreground zodiacal 
light, respectively. 
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Fig. 7. Top: The difference between the flux towards 
HH108MMS and the adjacent background in a cut along 
the major axis. The major axis goes through the position 
of maximum absorption in HH108MMS, where the offset 
is zero and T2 = T2.max, along a line connecting to the 
center of Qf . Bottom: The optical depth is derived after 
Eq.(|. 



observing the background than on the source. The coad- 
ded images at each raster position were projected on the 
sky and corrected for field distortion to derive the final 
mosaics. The mosaics have a pixel scale of 3". The dif- 
ferent raster maps of the same polariser are registered on 
the brightest pixel of IRAS 1833 1-0035. We measure a rel- 
ative shift of the source of < 1 pixel from one cycle to 
the next. The polarised signal is found to be consistent 
between the cycles, so that an average of all cycles gives 
the final mosaic image for each polarisers. 
From each mosaic we subtract the background. Its value is 
estimated on the same sky area. The Stokes elements are 
calculated following the notation and calibration factors 
given by Siebenmorgen (1999). The background at 12/im 
and 14/im is mostly due to zodiacal light. By assuming the 
latter to be unpolarised, we get an internal calibration of 
the instrumental polarisation which is consistent with the 
numbers cited in Siebenmorgen (1999). 
Fig. I displays the 14^m fluxes of IRAS18331-0035 as 
measured through the three polarisers PI, P2 and P3 over 
three observing cycles. For a given polariser, the fluxes 
agree within the uncertainty limits. However, the differ- 
ent polarisers have different fluxes, so the mid infrared 
emission of IRAS18331-0035 is polarised. Fig. | is simflar 
and refers to the 12/im fllter. Altogether we find: 

- 12/im: p = 28.9 ± 4.6%, 6i = 115 ± 5°. 

- 14/im: p = 15.2 ± 4.8%, 6i = 108 ± 9°. 

Repeating the observing procedure on either an unpo- 
larised standard star (HIC085371) or the zodiacal light. 



the fluxes through all polarisers are implying that the in- 
strumental polarisation is smaller than 1.5% (Siebenmor- 
gen 1999). 

3. Results 

3. 1 The optical depth profiles of the protostars 

To derive the density structure of the protostellar systems, 
HH108MMS and Ql, we transform the mid infrared sur- 
face brightness into an optical depth. The underlying ge- 
ometrical configuration is depicted in Fig. |^. HH108MMS 
is embedded in the Serpens molecular cloud and viewed 
from Earth through a veil of zodiacal light. We distinguish 
four regions that contribute to the observed flux: 

— The background with radiation of intensity Iq that 
heats the molecular cloud from behind. 

— The far and near side of the molecular cloud. They 
have constant optical depths ri , T3 and constant tem- 
peratures Ti,T3, respectively. 

— The absorbing core itself with corresponding parame- 
ters T2 and T2. 

— The interplanetary dust which is tenuous and emits 
zodiacal light of intensity I^. 

To simplify the analysis, we assume: 

— Ti — T3, the absorbing core is located in the optical 
center of the molecular cloud; 
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— Ti = T3 , the temperature is uniform over the molecular 
cloud; 

— Bi4f^n'L{Ti) = 0, the molecular cloud is cold and has no 
PAH emission; 

— -Bi4/im(?2) = 0, the absorbing core is isothermal and 
has also no mid infrared emission of its own. 

In the direction of the absorbing core (source) we can thus 
write 



/s = + /n • e-(2-i+-2) 



(1) 



Looking at a position next to the absorbing core, one re- 
ceives the background intensity 



h = h + he 



(2) 



If the zodiacal light is constant on scales of a few arcsec, 
this leads to 



h - Ih = In ■ e 



-2ri 



(3) 



Of course, the optical depth T2 varies across the absorbing 
core. Let us denote its and by 

-fs.max the observed intensity towards that position. An es- 
timate for T2.max IS obtained from the 1.3mm peak flux in 
a 11" beam of S'l.smm = 175 mJy (Chini et al. 1997) and 
the dust temperature of 13 K (Chini et al. 2000). From the 
formula = Bi,{T) fl with a solid angle = 2.2 • 10^^, 
we find that Ti.amm ~ 7 • 10~^. This is an average over the 
11" beam. As the source is centrally peaked at 1.3mm, we 
expect the optical depth in a 2" x 2" area, correspond- 
ing to the pixel size of the ISOCAM map, to be a few 
times larger. When we adopt a dust extinction coefficient 
of Ki.3mm = 0.02 cm^ per g of interstellar matter, as would 
be appropriate for a class object (Kriigel & Siebenmor- 
gen 1994), and the total 1.3mm fiux of 282mJy, we derive 
a total gas mass of O.SM©. 

To convert ri amm into the maximum optical depth at 
14/zm, T2,max, wc employ the dust models by Kriigel & 
Siebenmorgen (1994). For the diffuse interstellar medium, 
we read off from their Fig. 12 an optical depth ratio 
■''i4;im/'''i.3mm 1900. For grains of the same size in cold 
and dense cores where they have ice mantles and a fluffy 
structure, the ratio is about 1000. We thus end up with an 
estimate for T2,max that is above 10 and thus comfortably 
above unity. We may now safely assume that cxp (— T2,max) 
is small so that (see Eq.(3) 



h 



(4) 



Inserting Iqc^^'^^ into Eq.(||) gives the following approxi- 
mate expression for the variation of the optical depth in 
a cut through the major axis of the core. 



T2 



In 1 



Ih - Is, 



(5) 



It is independent of the optical depth of the molecular 
cloud, 2ti. The background intensity /b is determined as 
an average over two strips that run parallel to the major 
axis, one above, the other below it; we find /b = 4.3 • 
10~^^erg/s cm^ Hz ster. 

We can also derive Iq from Eq.(^ and thus the rate at 
which the cloud is heated from outside if we can guess ti . 
A rough number for 2ri in Eq.(^) comes from the H2CO 
line observations in the Serpens cloud by Ungerechts & 
Glisten (1984) who propose a visual extinction of 16 mag. 
With a ratio of rv/Ti4^m ~ 50 for dust in the diffuse 
medium, one obtains 2ti ~ 0.3. Our ISOCAM data then 
yield an intensity at 14/im: Iq ~ e'^-^ ■ 30/i Jy arcsec^^ or 

/o(14^ni) = 1.7 • 10~^'^erg/s cm^ Hz stcr, 

which is about two times higher than in the solar vicinity 
where /isRF(14/^m) = 7- 10^^^ erg/s cm^ Hz ster (Perault 
1987). So the value of Iq is reasonable considering that it 
refers to a star forming region where the radiation field is 
stronger. 

The combination of the ISOCAM and 1.3mm images al- 
lows to estimate the extinction ratio Ki4;im/Ki.3mm- When 
we roughly evaluate from the ISOCAM image the average 
optical depth of the absorbing core over the 11 beam of 
the 1.3mm map, we find for this average 3 (Fig. 
Therefore, 



^14^1111 _ (t2 + 2ti) 



14^m 



470 



Kl.3r 



(6) 



3.2 The density structure of the protostars 



From the variation of the optical depth over the source 
(Fig. 0) one can derive the density structure of the proto- 
stellar core. Consider a spherical cloud of Radius R and 
density p{r) — pQr~°'. The optical depth T2{x) through 
the cloud at an offset x < R from its center is given by 



T2 (x) = K po 



(7) 



where k is the dust absorption coefficient. For determining 
the exponent a of the density distribution, the product 
K/9o is an irrelevant proportionality factor. The value of 
T2 (x) depends on the upper boundary of the integral and 
thus on R. From the upper panel of Fig. 0, we estimate a 
source size for HH108MMS of 60" which corresponds to a 
cloud radius R = 10000 AU. 

The variation of the optical depth and the best fit den- 
sity distributions are shown for the major and minor 
axis of HH108MMS in Fig. | and Fig. |, respectively. In 
both cuts we derive good fits for a density distribution 
p{r) oc 7--i-8±o.i^ fpj^g source Ql is more difficult to ana- 
lyze because of the vicinity of HH108MMS and the IRAS 
source. 
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major axis 



radius are 10^^ and 10^^ cm, there is a visual extinction to 
the center of 220 mag, the total mass and total luminosity 
are 1.1 and 2.5L0, respectively. The spectral energy 
distribution is shown in Fig. 
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Fig. 9. The optical depth profile along the minor axis of 
HH108MMS (symbols). The solid lines show the variation 
of optical depth in a spherical cloud model with a power 
law density distribution, p oc 
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3.3 The spectral energy distribution of IRAS18331-0035 

We do a simple spherical radiative transfer calculation 
for IRAS18331-0035 to derive estimates for the source 
structure; the dust model is from Kriigel & Siebenmor- 
gen (1994) for fluffy and compact particles, the numerical 
code is described in Siebenmorgen et al. (1992). The few 
photometric points between 12//m and 1.3mm are fit by a 
source with a l/r~density distribution. Its inner and outer 
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Fig. 8. The optical depth profile along the major axis 
of HH108MMS as from Fig. | (symbols). The solid lines 
show the variation of optical depth in a spherical cloud 
model with a power law density distribution, p oc r^". 



3.4 Mid-infrared polarisation images 

Polarisation images are presented in Fig. |4| and Fig. |[ 
They have the following characteristics: 

- Mid-infrared polarisation is detected, but only along 
the absorbing filament. 

- The polarisation is strongly associated with the sources 
HH108MMS, Ql and IRAS18331-0035. 

- The degree of polarisation is high and reaches 50%. 

- Except for the western part of IRAS18331~0035, there 
is definitely some alignment between the polarisation 
vectors and the filament. 

- Because of their low temperatures, HH108MMS and 
Ql do not emit at 14/im, so polarised emission can be 
ruled out. 

- Because the background light source, as viewed from 
us, is behind the protostars, we exclude the possibility 
of polarisation from scattering. Furthermore, scatter- 
ing would require micron-sized dust particles. 

3.5 The degree of polarisation in the mid-infrared 

We now estimate whether an extinction optical depth at 
lAivcsi of order two can polarise the light by 30% or more. 
This has to be checked because polarisation through ex- 
tinction in the visible is always much smaller. In the visi- 
ble, we have the empirical relation pv ^ 0.03 ry between 
the visual optical depth ry and the degree of polarisation 
pv- 

We therefore compute the cross section of spheroids in the 
Rayleigh limit when the wavelength is much bigger than 
the particle. A spheroid has three major axes, a > b > c. 
If it has the shape of a cigar (prolate spheroid), a > b = c; 
if it is like a pancake (oblate spheroid), a = b > c. The 
absorption cross section at wavelength A of a spheroid of 
volume V when the electric field vector is parallel to axis 
a is 



C„ = 



2ttV 



Im 



1 



1 



(8) 



Im{z} signifies the imaginary part of a complex number 
z] e(A) denotes the complex dielectric permeability; La 
is a shape factor which lies between zero and one and 
depends on the eccentricity of the grain (see, for example, 
Bohren & Huffman 1983). There are analogous definitions 
for Cfc and Cc- For cigars. La + 2Lc = 1, for pancakes, 
2La+L,, = 1. 
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Fig. 10. The spectral energy distribution of IRAS18331- 
0035 

Suppose the spheroids are ahgned and rotate about the 
axis of greatest moment of inertia, which is axis & or c for 
cigars, and axis c for pancakes. Let the hght propagate 
in a direction perpendicular to the rotation axis and let 
C'e II rot:C'EJ_rot dcnotc the time-averaged cross sections 
in the case when the electric vector of a linearly polarised 
incident wave is parallel and perpendicular to the rotation 
axis, respectively. 

For spinning cigars, the cross section changes periodically. 
Because the mean of cos^ x over a rotation cycle is one 
half, one finds 

CE±rot = I [Ca + Cc] , Ce || rot = Cfe • (9) 

If a pancake rotates about its axis of maximum moment 
of inertia, which is c, no averaging is necessary and 

C'e _L rot = Ca , Ce II rot = Cc . (10) 

The optical depth r of a cloud for linearly polarised light 
changes with the direction of the electric field vector. Let 
■''max be the maximum value of r and Tmin its minimum. 
The corresponding electric field vectors are perpendicular 
to each other. For both cigars and pancakes, Tmax is pro- 
portional to CEirot and Taun to Ce II rot with the same 
proportionality factor. 

When unpolarised background light of intensity /* tra- 
verses the cloud, it is weakened to the observed intensity 




wavelength [/.im] 



Fig. 11. The percentage of polarisation due to perfectly 
aligned spinning spheroids of silicate and amorphous car- 
bon of axial ratio c/a =1/3 and c/a — 2/3. The effective 
optical depth r^g has the same value at all wavelengths. 
Solid lines refer to cigars, dotted to pancakes; in each box, 
the lower pair of curves to r = 1, the upper to t = 3. The 
size of the grain is irrelevant as long as it is much smaller 
than the wavelength. 

Putting Jobs = I* e^^""' defines an effective optical depth 

Tcff = - In > . (12) 

The degree of polarisation is obviously defined by 

P(A) = ^ , ■ (13) 

In case of weak extinction, Toff — ^ (T,„ax + Tmin) and p = 

5 (Tmax - Tnrin). If Tmax - Tmin is large, the polarisation 
goes to unity. 

In Fig. |ll| we calculate p{X) assuming an effective optical 
depth TeflF at all wavelengths of either 1 or 3. The frequency 
dependent dielectric permeability for silicate is from Laor 

6 Draine (1993) and for carbon from Zubko et al. (1996, 
their type BE) . We see that with an effective optical depth 
of two, even particles of very moderate elongation (a/c = 
3/2) produce very substantial polarisation, provided they 
are well aligned. Therefore, the observed high degree of 
polarisation at 14/im does not seem to pose a principal 
problem. 
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4. Discussion: Pre stellar evolution 

Among many other models, the self-similarity solutions 
of Shu et al. (1977) for the early evolution of protostars 
predict in the absence of rotation and magnetic fields a 
p oc r"^ profile in the outer envelope and a less steep dis- 
tribution p cx r~^l'^ in the inner part. Basu & Mouschovias 
(1995), who include magnetic fields and rotation, also find 
a power law density distribution in the envelope with an 
exponent between —1.5 and —1.8 and a constant den- 
sity central region of ~50AU size. Masunaga & Inutsuka 
(1999) perform calculations for the early collapse (before 
dissociation of H2 in the core) of a 1 star. The lumi- 
nosity stays below 0.1 L©, there is a core of 5 AU and the 
density goes like in the envelope. Their models should 
be applicable to HH108MMS and Ql of the present paper 
which also have no mid-IR embedded source. 

In a few cases, the predicted theoretical density profiles 
could be checked by 1.3mm observations of the dust emis- 
sion (Ward-Thompson et al. 1994, Andre et al. 1996, 
Ward-Thompson et al. 1999) and were generally corrob- 
orated. In this method, the conversion of fluxes into ab- 
solute column densities depends on the dust absorption 
coefficient at 1.3 mm and the grain temperature, whereas 
the power law exponent of the density distribution is sen- 
sitive to the temperature gradient. 

On the other hand, our derivation of the optical depth and 
the ensuing density profile from absorption measurements 
is temperature independent. ISOCAM observations simi- 
lar to ours were also carried out by Bacmann et al. (2000). 
For the sources which they present, the maximum optical 
depth at the cloud center is much smaller than those for 
HHIOSMMS and Ql, so we probe deeper into the cloud, 
still the density profiles in the envelopes roughly agree 
with what we derive. In the cloud center (< 2000AU), 
they find a flattening of the column density which we do 
not see. Such a flattening is expected to occur in the cen- 
tral lOOOAU of an isothermal sphere with temperature 
T ~ lOK, radius of lOOOOAU and mass of IM© simply by 
solving the hydrostatic equation (Bonnor 1956). 

The measured mid infrared to millimeter dust extinction 
ratio of Ki4/im/Ki3oopm ~ 470 should be compared with 
dust models of protostellar environments. The ratio is a 
factor 4 lower than expected for dust in the diffuse medium 
and it indicates that the grains in the dense and cold en- 
vironment of the protostar HHIOSMMS are, as expected 
(Ossenkopf 1993), of rather fluffy and composite nature. 
The measured mid infrared to millimeter dust extinction 
ratio of HHIOSMMS is already tending towards somewhat 
elongated grain structures: In the fluffy composite but 
spherical dust model by Kriigel & Siebenmorgen (1994) 
a ratio of Ki4^m/Ki300/im ~ 1000 is found. Because elon- 
gated particles are much better antennas at 1.3mm, they 
can give a lower value. A proof that the dust in the ab- 
sorbing cores and the IRAS source is indeed of elongated 



structure is found by the ISOCAM polarisation measure- 
ments. 

As the filament is seen in absorption against the back- 
ground and assuming that the majority of the grains are 
still sub-micron sized particles, one may neglect dust scat- 
tered light at 14/im. 

Consequently the most plausible mechanism to pro- 
duce the polarisation is dichroism. For polarisation due 
to extinction of elongated spinning and aligned dust par- 
ticles, the polarisation vectors (Figs. ^ and ||) indicate 
the magnetic field direction. It appears that the magnetic 
field vectors are roughly aligned with the absorbing fila- 
ment. Such ordered fields are also detected from S50/xm 
maps of thermal dust emission of prestellar cores (Ward- 
Thompson et al. 2000) and FIR polarimetry on the pro- 
tostar IRAS20503+6006 (Clements et al. 1999). 

5. Conclusion 

We study the protostellar system HHIOSMMS in the mid 
IR with ISOCAM. The 14.3/im (lw3) image reveals an 
extended (200" x 20") structure which is seen in absorption 
against the diffuse background radiation. Within it, we 
detect IRAS1S331-0035 emitting a flux of a few mJy and 
at least two absorbing cores, of which one coincides with 
the protostar HHIOSMMS. 

We transform the 14/im surface brightness of the absorb- 
ing cores into optical depth profiles from which we derive 
the density structure of the protostars. Along the minor 
and major axis of the absorbing cores the density profile 
is consistent with a r~^'^ power law. 
From the comparison to the 1.3 mm map (Chini et 
al. 1997) we derive observationally an extinction ratio 
/Ki3oo/jm ~ 470. When compared to dust models 
(Kriigel & Siebenmorgen 1994), this ratio points towards 
rather fluffy and elongated grain structures, as expected 
for cold and dense environments. 

The structure of IRAS1S331-0035 is estimated by apply- 
ing a radiative transfer model to fit its spectral energy 
distribution. A reasonable fit is obtained using a l/r den- 
sity structure, a total cloud mass of Afgas = 1. IM© and 
a luminosity of L = 2.5 Lq. The model envelope has a 
visual extinction of Ay ~ 220mag, where we measure Ay 
from the cloud surface to the inner loci where grains are 
evaporating. 

We also present the first mid infrared polarisation images 
of protostellar absorbing cores prior to the formation on 
a luminous shock front. 

The polarisation pattern is well correlated with the triplet 
system and follows roughly the absorbing filament. The 
most plausible explanation for the polarisation is extinc- 
tion of the background radiation by elongated spinning 
and aligned dust particles. Therefore the polarisation vec- 
tors indicate the magnetic field structure of the protostel- 
lar triplet system. 
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We detect a fractional polarisation which is much higher 
(>15%) than known for optical wavelengths (< 10%). For 
various kinds of rotating dust particles it is shown that 
dichroic polarisation can indeed produce such high degree 
of mid infrared polarisation. 
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